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^ ' Abstract. A new approach to the study of the transition point in a class of two dimensional Wess- 

[ Zumino models is presented. The method is based on the calculation of rigorous lower bounds on the 

. ground state energy density in the infinite lattice limit. Such bounds are useful in the discussion of 

' supersymmetry phase transition. The transition point is then determined and compared with recent 

, results based on large-scale Green Function Monte Carlo simulations with good agreement. 

^ ■ 

o ■ 

I The simplest theoretical laboratory to study non perturbative dynamical supersym- 

■ metry breaking is the two dimensional = 1 Wess-Zumino model that involves chiral 

' ! superfields with no vector multiplets. Let us remind the (continuum) N = I supersym- 

Q I metry algebra, {Qa, Qp} = 2(/'C)c^j8. Since Pj are not conserved on the lattice, a lattice 

^ ■ formulation of a supersymmetric model must break this algebra explicitly. A very im- 

> ! portant advantage of the Hamiltonian formulation is the possibility to conserve exactly 

k3( I a key subalgebra of this relation; specializing to 1 + 1 dimensions, in a Majorana basis, 

5-2 ; 7o = C = C72, 71 = i^3, the algebra becomes Qj = Ql=pO = H, {21,22} = = 

. . 2P. On the lattice, since H is conserved but P is not, we can pick up one of the super- 
charges, say, QI = H, build a discretized version Qi and define the lattice Hamiltonian 

H = Ql. 

The explicit lattice model is built by considering a spatial lattice with L sites. On 
each site we place a real scalar field % together with its conjugate momentum 
such that [pn, <p,n] = —idn.m- The associated fermion is a Majorana fermion x^fa.n with 
a = 1,2 and {\ira,n,Wb,m} = 8a,b^n,m, v4,« = V^","- The discretized supercharge Ql = 
L^^i [Pnm,n - ( ''^"+'2''^"'' + V^2,«], with arbitrary V{(p) (called prepotential) can 

be used to define a semipositive definite lattice Hamiltonian H = Q\. Notice that Q\ = H 
is enough to guarantee that Eq = (^ol^l^o) > 0, that all eigenstates of H with E > 
are paired in doublets and that £0 = if and only if supersymmetry is unbroken, i.e., the 
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ground state is anihilated by Qi. 

Rigorous results from the continuum can be found in [1]. On the lattice, accurate nu- 
merical results are available [2, 3], although a clean determination of the supersymmetry 
breaking transition remains rather elusive. All predictions and results indicate: for the 
model with cubic prepotential, V = (p^, unbroken supersymmetry; for the model with a 
quadratic prepotential, V — ^2^^ + ^, dynamical supersymmetry breaking. Along a line 
of constant A2 the results for a quadratic potential indicate the existence of two phases: a 
phase of broken supersymmetry with unbroken discrete Z2 at high Aq and a phase of un- 
broken supersymmetry with broken Z2 at low Ao, separated by a single phase transition. 
On the other hand, strong coupling expansion demonstrate that for a polynomial V{(p), 
the relevant parameter is just its degree q. For odd q, strong coupling expansion and 
tree-level results agree and supersymmetry is expected to be unbroken. This conclusion 
gains further support from the nonvanishing value of the Witten index [4]. For even q in 
strong coupling expansion, the ground state has a positive energy density also for L ^ 00 
and supersymmetry appears to be broken for all Aq. On the other hand, weak coupling 
expansion predicts unbroken supersymmetry when Aq < [2]. 

We used two different approaches to investigate the pattern of dynamical supersym- 
metry breaking in a class of two dimensional Wess-Zumino models. In the first one, 
[2, 3], the numerical simulations were performed using the Green Function Monte Carlo 
(GFMC) algorithm and strong coupling expansion. The GFMC is a method that com- 
putes a numerical representation of the ground state energy density on a finite lattice 
with L sites in terms of the states carried by an ensemble of K walkers. By perform- 
ing numerical simulations along a line of constant A2 = 0.5, we determined the nu- 
merical value of Ao separating a phase of broken supersymmetry from a phase of un- 
broken supersymmetry. The usual technique for the study of a phase transition is the 
crossing method applied to the Binder cumulant [2]. The crossing method consists in 
plotting B vs. Ao for several values of L. The crossing point Xq{Li,L2), determined 

by the condition B{Xq,Li) — B{Xq,L2) is an estimator of Aq'^^. The value obtained is 

Aq ' = —0.48 ± 0.01. The main source of systematic errors in this method is the need to 
extrapolate to infinity both K and L. For this reason, an indepentend method to test the 
numerical results of [2] is welcome. 

The second method is based on a new approach to the study of the supersymmetry 
phase diagram introduced in Ref. [5] and is based on the calculation of rigorous lower 
bounds on the ground state energy density in the infinite lattice limit. The Hamiltonian 
lattice version of the Wess-Zumino model conserves enough supersymmetry to prove 
that the ground state has a non negative energy density p > 0, as its continuum limit. 
Moreover, the ground state is supersymmetric if and only if p = 0, whereas it breaks 
(dynamically) supersymmetry if p > 0. Therefore, if an exact positive lower bound Plb 
is found with < Plb < p, we can claim that supersymmetry is broken. 

The relevant quantity for our analysis is the ground state energy density p evaluated on 

the infinite lattice limit, p — limi^^ — limi^oo p{L). It can be used to tell between 
the two phases of the model: supersymmetric with p = or broken with p > 0. We have 
shown in [5] how to build a sequence of bounds p^^^ which are the ground state energy 
density of the Hamiltonian H with modified couplings on a cluster of L sites. We now 




explain how to exploit the sequence of bounds p^^^ to determine the phase at a particular 
point in the coupling constant space. In order to do so, we compute numerically p*^^) at 
various values of the cluster size L. Ifp^'^) >Ofor some L, we can immediately conclude 
that we are in the broken phase. On the other hand, if we find a negative bound we cannot 
conclude in which phase we are. However, we know that p^^) — > p for L — > oo and the 
study of p as a function of both L and the coupling constants permit the identification 
of the phase in all cases. To test the method we studied in detail the case of the quadratic 
prepotential and discussed the dependence of p on Aq at a fixed value of X2 — 0.5 (as in 
Ref. [2]). 

On the left side of Fig. 1 a reasonable qualitative pattern of the curves representing 
P^^\Aq) is shown. Notice that a single zero is expected in p(^)(Ao) at some Ao = 
Xq{L). Since lim/^^cop*^^-* = p, we expect that Ao(L) — > Xq for L — > 00, allowing for 
a determination of the critical coupling Xq. The results of the energy for the lower 
bound p^^) (Ao) for all cluster sizes behaves as expected: it is positive around Ao = and 
decreases as Ao moves to the left. At a certain unique point Aq (L), the bound vanishes 
and remains negative for Aq < Aq (L). This means that supersymmetry breaking can be 
excluded for Ao > min^Ao (L). Also, consistency of the bound means that Xq{L) must 
converge to the infinite-volume critical point as L — > 00. Since the difference between the 
exact Hamiltonian and the one used to derive the bound is ^( 1 /L) , we can fit Aq (L) with 
a polynomial in 1/L. This is shown on the right side of Fig. 1 [5], where also the GFMC 
result is quoted. The best fit with a parabolic function gives Aq = —0.49 ± 0.06, quite in 
agreement with the previous Aq Qpj^^ = —0.48 ± 0.01 using the GFMC algorithm. 
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